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LONG-TERM  GOALS 

Sound  propagation  in  the  areas  of  continental  shelves  is  complicated  due  to  many  three-dimensional 
(3-D)  oceanographic  and  marine  geologic  features,  such  as  shelfbreak  frontal  systems,  nonlinear 
internal  gravity  waves  and  topographic  variability.  The  long-term  goals  of  this  project  are  targeted  on 
understanding  the  3-D  sound  propagation  effects  caused  by  these  environmental  factors,  and  also  on 
applying  the  3-D  sound  propagation  physics  to  acoustic  inversions. 

OBJECTIVES 

Due  to  the  horizontal  inhomogeneity  in  the  water  column  caused  by  shelfbreak  frontal  systems, 
nonlinear  internal  gravity  waves  and  other  physical  oceanographic  processes,  horizontal  refraction  of 
sound  occurs  and  produces  significant  3-D  acoustic  propagation  conditions.  One  of  the  objectives  of 
this  project  is  to  develop  an  efficient  and  accurate  3-D  sound  propagation  model  (both  theoretical  and 
numerical  models)  for  applications  related  to  low-frequency  sonar  systems. 

In  addition  to  the  3-D  sound  propagation  study,  this  project  also  has  the  investigation  on  the  feasibility 
of  two  specific  inversion  techniques  in  the  area  of  continental  slopes:  (a)  source  localization  and  (b) 
bottom  geoacoustic  inversions.  This  work  requires  a  solid  understanding  on  3-D  sound  propagation, 
and  it  is  closely  connected  to  the  first  topic. 

APPROACH 

The  technical  approaches  employed  in  the  first  study  (3-D  sound  propagation)  include  theoretical 
analysis,  numerical  computation  and  real  data  analysis.  A  3-D  normal  mode  method  is  used  to  study 
canonical  environmental  models  of  shelfbreak  front  systems  and  nonlinear  internal  wave  ducts.  The 
WHOI  3D  Parabolic-Equation  (PE)  wave  propagation  model  [1]  and  a  vertical-mode  horizontal-PE 
model  of  the  Weinberg-Burridge  formalism  [2]  are  used  to  study  sound  propagation  in  real  situations. 
The  idea  of  normal  mode  back-propagation  is  implemented  for  the  acoustic  inversion  work.  A  two- 
dimensional  (2-D)  approach  is  developed  first,  and  then  will  be  generalized  for  3-D  cases. 
Experimental  data  collected  from  recent  work  on  the  New  Jersey  shelf  (SW06  [3]),  the  East  China  Sea 
(QPE  [4])  and  the  South  China  Sea  (ASIAEX  [5]  and  NLIWI  [6])  are  analyzed  with  collaboration  with 
experiment  participants. 
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WORK  COMPLETED 


The  tasks  completed  in  the  year  are  described  below. 

1.  3D  PE  model  development 

The  WHOI  3D  PE  program  has  been  modifed  for  parallel  CPU  and  GPU  computing.  Also,  a.  grid  size 
requirement  is  also  derived  using  the  sampling  theory: 


(  K  max  =  n/Ay,Aky  =  n/Ly  (|y|  <Ly,\ky\<  ky  max  ) 
DFT sampling  theory  j  kzmax  =  n/Az,Akz  =  n/Lz  (|z|  <  Lz,\kz\  <  kzmax) 

Uemax  =n/A6,Ak0  =n/Lg  (\d\  <  Le,\ke  \  <  ke  max) 


where  ky>Zte >are  the  wavenumbers  along  the  transverse,  vertical  and  angular  axes.  This  relations  are 
fundamental,  and  they  actually  govern  the  requirements  for  the  model  grid  size.  Consieder  that  the  goal 
of  sound  propagation  modeling  is  to  resolve  all  of  the  arrivls  less  than  a  given  arrival  angle  \Q  <  fmax 
with  a  requried  angular  resolution  Afmin.  From  fmax  we  can  find  that  the  maximal  wave -numbers  need 
to  be  kyjZmax  >  kosin£max  and  kdmax  >  k0rsmfmax,  or  equivalently  for  Ay,  Az,  and  AP. 


Ay,Az  <  (sin^J  1  l0/2  and  Ad  <  (r  sin^J  1 (2  . 


On  the  other  hand,  from  the  requiement  of  angular  resolution  Afmin,  we  can  find  that  the  wavenumber 
increment  Aky  z  <  kocos£  A£min  and  Ako  <  ko  (rcosC A£mi„  +  sin<f  Ar),  or  equivalently  for  the  apertures  Ly, 
Lz  and  Lq, 


LyiZ  >  (cos<A<min)  1k0/2  and  Le  >(rcos<A<min  +sin<Ar)  1  k0/2  . 


From  this  analysis,  one  can  see  that  in  cylindrical  PE  the  models  resoultion  is  degrading  as  the  range 
goes  further.  To  overcome  the  problem,  two  improved  model  grids  are  suggested,  as  shown  in  Figure 

1.  The  performance  of  different  PE  models  is  shown  in  Figure  2,  and  one  can  see  that  the  Cartesian  PE 
and  the  cylindrical  PE  with  adaptive  grids  agree  with  each  other  very  well  on  the  solution  plane.  On 
the  other  hand,  the  coarse-grid  cylindrical  PE  produces  large  errors  due  to  the  degration  of  the  model 
resolution. 

2.  3-D  sound  propagation  effects  caused  by  topographic  variability 

Two  types  of  topographic  features  that  will  cause  3-D  sound  propagation  effects  are  studied  in  the 
year.  The  first  feature  is  the  seafloor  scours.  Because  the  bottom  scours  often  have  strong  directivity, 
propagating  sound  will  be  focused  along  the  scouring  direction  due  to  the  transverse  gradient  on  the 
bottom  depth.  Also,  the  horizontal  focusing  will  vary  with  the  source  position.  Using  the  method  of 
verical  model  and  horizontal  PE,  a  3-D  sound  propagation  model  with  the  adabaitc  mode  assumption  is 
built.  This  model  has  been  applied  for  an  SW06  example.  The  model  clearly  shows  that  due  to  the 
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seafloor  scours  in  the  New  Jersey  shelf  area  propagating  sound  can  be  focused,  as  shown  in  the  upper 
panel  of  Figure  3.  Two  source  positions  are  modeled  in  the  example,  and  the  resultant  focusing 
patterns  are  slightly  different.  Although  the  effects  of  the  seafloor  scours  are  profound,  they  will  be 
masked  episodically  by  the  presence  of  water  column  fluctuations,  such  as  nonlinear  internal  waves. 
Twenty-two-days  data  are  used,  and  a  long  time  average  is  taken  for  removing  the  episodic 
fluctuations.  The  data  shows  about  5-dB  intensification  due  to  the  horizontal  focusing.  The  data-model 
comparison  is  shown  in  the  lower  panel  of  Figure  3,  and  has  very  good  agreement. 

The  second  bathymetric  feature  that  will  cause  strong  3-D  sound  focusing  is  the  submarine  canyon.  A 
realistic  model  has  been  implemented  to  model  3-D  sound  propagation  in  the  North  Mien-Hua  canyon 
northeast  of  Taiwan,  see  Figure  4.  Because  of  the  concave  bathymetry,  the  sound  reflected  the  seafloor 
will  be  focused,  as  shown  in  the  lower  panel  of  Figure  4.  This  also  provides  a  plaussible  explanation 
for  the  intensified  ship  noise  recorded  at  a  hydrophone  array  during  the  QPE  experiment. 

3.  Whispering  gallery  modes  in  nonlinear  internal  wave  ducts 

Using  3-D  normal  mode  theory,  general  solutions  of  sound  pressure  field  in  a  curved  internal  wave 
duct  is  found: 


p(r,  6,  z)  = 


1' 


Ad  e*(fc2(z)~£m)2z  +  au  e-*(fe2(z)-£m)2z 

vertical  modes 


X  [Bout  «mr)  +  Bin  H™„  «mr)]  X 

radial  modes 


where 


=/„(*) +  ;r„  (*) 

H^\x)=Ux)-iYv(x) 

The  field  is  decomposed  into  vertical  and  radial  components.  For  each  vertical  mode,  there  is  a  set  of 
radial  modes,  which  is  a  form  of  Bessel  functions  desorbing  the  horizontal  component  of  the  field. 
Also,  the  radial  modes  can  propagate  along  the  wave  front.  Two  types  of  modes  are  found.  One  is  the 
whispering  gallery  modes  formed  by  the  sound  energy  trapped  by  the  outter  wave.  The  other  is  full 
bouncing  mode  generated  by  the  sound  energy  bouncing  between  inner  and  outter  waves. 

An  exmple  of  50  Hz  sound  propagation  in  a  curved  nonlinear  internal  wave  duct  with  25  km  curvature 
is  shown  in  Figrue  5.  The  upper  panel  shows  the  vertical  modal  structure,  and  the  lower  panels  are  the 
radial  modes  associated  with  the  second  vertical  mode.  From  the  radial  mode  shape,  one  can  see  that 
the  radial  mode  1  is  a  whiserping  gallery  mode  since  it  is  formed  by  the  sound  energy  trapped  by  the 
outter  wave,  resulting  a  single  hump  close  to  the  outter  boundary. 

4.  Adapative  normal  mode  back-propagation  approach 

A  variety  of  localization  methods  with  normal  mode  theory  have  been  established  for  localizing  low 
frequency,  broadband  signals  in  a  shallow  water  environment.  Gauss-Markov  inverse  theory  is 
employed  in  this  paper  to  derive  an  adaptive  back-propagation  approach.  Joining  with  the  maximum  a 
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posteriori  mode  filter,  this  approach  is  capable  separating  signals  from  noisy  data  and  back- 
propagating  them  without  significant  influence  from  noise.  Numerical  simulations  are  presented  to 
demonstrate  the  robustness  and  accuracy  of  the  approach  presented,  along  with  comparisons  to  other 
methods.  Applications  to  real  data  collected  from  the  SW06  experiment  are  presented  in  Figure  6,  and 
the  effects  of  water  column  fluctuations  with  scales  from  nonlinear  internal  waves  to  shelfbreak  front 
variability  are  observed. 

RESULTS 

An  improved  3-D  cylindrical  PE  model  with  consistent  azimuthal  resolution  is  developed.  This 
cylindrical  model  can  complement  the  3-D  Cartesian  PE  for  a  cylindrical- wave  like  field.  (The  3-D 
Cartesian  PE  is  not  good  for  a  cylindrical-wave  like  field  because  its  PE  approximation  errors  increase 
in  greater  azimuth. 

Reaslitic  models  are  implemented  to  study  the  bathymetric  effects  on  3-D  sound  refraction  and 
focusing.  Experimental  data  also  demostrate  significant,  nonnegotiable  3-D  effect.  Specifically,  the 
cancave  seafloor  of  a  submarine  canyon  can  cause  sound  focusing,  more  than  10-dB  intensificiation  is 
expected.  Study  on  the  whispering  gallery  effect  in  curved  internal  wave  ducts  indeed  provides  a 
simple,  but  elegant,  theoretical  model  to  the  problem. 

Finally,  an  adaptive  normal  mode  back-propagation  approach  for  low-frequency  broadband  sound 
source  localization  in  a  shallow-water  ocean  is  established.  Gauss-Markov  inverse  theory  is  used  in 
both  mode  filtering  and  back-propagation.  This  unifies  the  adaptive  normal  mode  back-propagation 
approach.  The  method  can  smoothly  adapt  to  the  signal-to-noise  ratio  and  ensure  an  optimal  balance 
between  robustness  and  accuracy. 

IMPACT/APPLICATIONS 

The  potential  relevance  of  this  work  to  the  Navy  is  on  increasing  the  capability  of  Naval  sonar  systems 
in  shallow  water  areas.  The  contributions  of  the  effort  on  studying  3-D  sound  propagation  effects  will 
be  on  assessing  the  environment-induced  acoustic  impacts.  In  addition,  the  investigation  of  acoustic 
inversions  directly  relates  to  the  Navy  sonar  operation. 

TRANSITIONS 

The  3-D  sound  propagation  model  has  a  potential  for  transition. 

RELATED  PROJECTS 

Experimental  data  were  collected  from  the  ONR  ASIAEX,  SW06  and  QPE  projects.  Also,  collabration 
with  Dr.  Megan  Ballard  of  ARL,  UT  at  Austin  is  established  for  the  study  on  the  horizontal  refraction 
of  propagating  sound  due  to  seafloor  scours. 
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(a)  Spatial  domain 


(b)  Wavenumber  domain 


Upper  panels:  Upsampling  of  angular  grid  points  in  the  cylindrical  PE  model 


Lower  panel:  A  fixed  arc-length  grid  for  the  cylindrical  PE  model 


Figure  1.  Two  improved  model  grids  for  the  cylindrical  PE  model 

[  Upper  panels  show  the  zero-padding  on  the  wavenumber  spectrum  results  in  upsampling  of 
angular  grids.  Lower  panel  shows  the  free  propagation  path  occurred  in  a  fixed  arc-length  grid.] 
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Modeling  comparisons 

Propagate  over  seamount,  off  center 
Source  at  250  m,  100Hz 
4  cases  -  (1)  Nx2D,  (2)  Cartesian,  (3)  coarse-grid 
cylindrical  and  (4)  adaptive-grid  cylindrical 


Nx2D  PE  solution,  frequency  100.00  Hz  dB 


Cartesian  3D  PE  solution,  frequency  100.00  HaB 
o 

(2) 


X<m) 

adaptive  azimuth  grid 

Cylindrical  3D  PE  solution,  frequency  100.00  Hz  c 


Upper  panels:  3D  PE  results 


Lower  panels:  3D  PE  differences.  The  baseline  is  the 
adaptive-grid  cylindrical  PE. 


Figure  2.  PE  model  comparisons 

[Sound propagation  over  a  seamount  are  computed  by  different  3D  PE  models,  including  (1)  Nx2D 
(2)  Cartesian,  (3)  coarse  cylindrical  grid  and  (4)  adaptive  cylindrical  grid.  The  computed  sound 
fields  are  shown  in  the  upper  panels,  and  the  model  differences  are  shown  in  the  lower  panels.  The 
Cartesian  PE  and  the  cylindrical  PE  with  an  adaptive  grid  have  very  good  agreement.  ] 
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(a)  NRL  Source  (b)  MSM  Source 
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Upper  panels:  the  horizontal  refraction  of  propagating  sound  from  two  different 
_ moored  sources  in  the  SW06  experiment.  (Model  results) _ 
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Lower  panels:  Data-model  comparison  of  the  horizontal  refraction  caused  by  the 
seafloor  scours  in  the  SW06  experiment.  The  data  is  shown  as  probability 
density,  and  the  model  is  shown  in  blue  curves  and  presenting  the  mean  field. 


Figure  3.  Horizontal  refraction  of  sound  caused  by  seafloor  scours 


[  Acoustic  data  and  models  demonstrate  the  acoustic  effect  of  seafloor  scours  on  horizontal 
refraction.  The  model  well  captures  the  energy  focus  seen  in  the  data. 

The  sound  frequencies  are  100 ,  200  and  300  Hz.  ] 
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(a)  Bathymetry  of  the  North  Mien-Hua  Canyon,  northeast  of  Taiwan 
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Depth-averaged  sound  intensity  over  top  20-m  water,  3D  model 
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Sound  intensity  structure  along  the  solid  line  in  panel  (c),  3D  model 
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(b)  Comparison  of  Nx2D  and  3D  PE  models.  Acoustic  frequency  300  Hz. 


Figure  4.  Modeling  of  sound  propagation  over  a  submarine  canyon 

[  The  bathymetry  of  the  canyon  system  is  shown  in  the  panel  (a),  and  the  PE  models  are  shown  in 
panel  (b).  From  the  differences  between  the  Nx2D  and  3D  PE’s,  one  can  see  that  the  canyon 
bathymetry  causes  very  strong  sound  focusing.  The  cylindrical  spreading  loss  is  removed  in  the  plots 

to  reduce  the  dynamic  range  of  the  TL  variability.  ] 
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(a)  Vertical  modes  of  50  Hz  sound  inside  and  outside  the  internal  waves 
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(b)  Radial  modes  of  50  Hz  sound  in  the  internal  wave  duct 


Figure  5.  Whispering  gallery  modes  in  a  curved  nonlinear  internal  wave  duct. 

[3-D  normal  mode  theory  is  used  to  study  the  whispering  gallery  modes  occurred  in  a  curved 
nonlinear  internal  wave  duct.  For  each  vertical  mode,  as  shown  in  the  panel  (a),  there  are  a  set  of 
radial  modes  describing  the  horizontal  component.  In  this  example,  an  internal  wave  duct  with  25 
km  radius  is  considered.  The  radial  modes  of  the  second  vertical  mode  are  shown  in  the  panel  (b). 

One  can  see  that  the  radial  mode  1  humps  near  the  outer  boundary, 
presenting  a  whispering  gallery  mode.  ] 
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MSM  source  localization  results 
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(a)  SW06  experimental  area  and  mooring  locations  (b)  Effects  of  nonlinear  internal  waves  on  MSM 

source  localization. 
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(c)  Distributions  of  MSM  source  localization  results  from  the 
normal-mode  back-propagation  approach.  The  true  source  location 
was  at  19.75  km  in  distance  and  55  m  in  depth. 


Figure  6.  Source  localization  using  the  adaptive  normal-mode  back-propagation  method 

[  The  adaptive  normal-mode  back-propagation  method  is  applied  to  the  SW06  data.  Processed 
signals  are  from  one  of  the  moored  sources  in  the  experiments,  the  MSM  source  shown  in  the  panel 
(a),  and  the  source  localization  results  are  shown  in  the  panel  (b)  to  correlate  with  the  nonlinear 
internal  wave  signals  measured  at  the  environmental  mooring  ENV#32.  A  positive  correlation  can 
be  seen.  Distributions  of  the  source  localization  results  are  shown  in  the  panel  (c),  showing  very 
good  source  range  estimation.  The  bimodal  distribution  of  the  source  depth  estimates  is  most  likely 

caused  by  unresolved  mode  coupling  effects.  ] 
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